The nonlinear damping of Bose-Einstein condensate oscillations at ultra-low 

temperatures 
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We analyze the damping of the transverse breathing mode in an elongated trap at ultralow 
temperatures. The damping occurs due to the parametric resonance entailing the energy transfer to 
the longitudinal degrees of freedom. It is found that the nonlinear coupling between the transverse 
and discrete longitudinal modes can result in an anomalous behavior of the damping as a function 
of time with the partially reversed pumping of the breathing mode. The picture revealed explains 
the results observed in fan]. 



The problem of the damping of oscillations in a Bose- 
condensate at T — and, at least, at ultra-low temper- 
atures within T < \i -C T c , \i being chemical potential 
and T c being Bose-condensation temperature, represents 
one of most interesting problems in the physics of Bose- 
Einstein condensation. The specific feature of the prob- 
lem is associated with an isolation of the system from 
the environment and, therefore, absence of any sort of 
thermal bath. In this aspect the problem has much in 
common with the problem of the loss of coherence and 
the origin of energy dissipation in an isolated quantum 
many-particle system. Vast majority of experimental |jl|]- 
H and theoretical efforts in this field are referred 

to the range of relatively high temperatures T 3> fo^i 
where u>i are the frequencies of a parabolic trap. Un- 
der these conditions normal excitations in essence form 
an interior thermostat. The damping in the case con- 
cerned is associated with the interaction of oscillations 
and thermal excitations accompanying by direct transfer 
of energy to the thermal cloud. 

In the work of authors [[l3) the problem is consid- 
ered for the damping of radial oscillations of the con- 
densate in an elongated trap at T = 0. The choice 
of the trap geometry was not casual. As is found in 
|0 for a 2D isotropic harmonic trap, the non-linear 
Schrodinger (Gross-Pitaevskii) equation has a scaling- 
invariant solution for the evolution of the condensate and 
non-condensate subsystems at an arbitrary variation of 
uj(t) in time. In this case the breathing mode in the 2D 
system does not decay irrespective of the magnitude of 
the amplitude. It was shown [jl5| that the unique prop- 
erties result from the presence of a hidden symmetry in 
the description of such 2D system. 

On the other hand, it is found in Jl3) that the trans- 
verse breathing mode can experience a specific damping 
in an elongated trap. The origin of such damping lies in 
the parametric resonance connecting the breathing mode 
with longitudinal sound excitations. The parametric res- 
onance arises from the radial oscillations of the density, 
resulting in oscillation of the longitudinal sound velocity. 
In fact, the results obtained give an appropriate descrip- 
tion only for the initial stage of damping, see below. 

Recently, it has been published very interesting ex- 



perimental investigation on the damping of the breath- 
ing mode in an elongated trap at ultra-low tempera- 
tures T < ji « T c p6| . The work contains a num- 
ber of remarkable results. For the limited initial am- 
plitude of oscillations characterized by the scale of rel- 
ative variation \SR±/R±\ ~ 0.03 — 0.04 for the trans- 
verse Thomas-Fermi radius, the authors have observed 
the slow and monotonous decay with the record quality 
factor Q = ^o/7o > 2000. Here ujq = 2cj^ and 7 are 
the frequency of breathing mode and the damping rate, 
respectively. The temporal scale of damping is about 1 s. 
However, provided the initial amplitude is increased by 
a factor of three, the whole picture changes drastically. 
Here on the temporal scale of 200 ms the oscillation am- 
plitude reduces by about three times. Next, the ampli- 
tude unexpectedly starts to grow, reaching the maximum 
of about one-half of the initial value, and then decays. 
The aim of the present study is to explain an anomalous 
picture of damping the breathing mode observed in [jl6l| . 

The standard theory of parametric resonance implies 
an invariability for the amplitude of oscillation of one 
of parameters in the system. In the general case this 
corresponds to the steady feeding of energy. In the iso- 
lated system the situation is different. The amplitude 
of breathing mode b± decreases in the course of the en- 
ergy transfer into the subsystem of longitudinal excita- 
tions due to parametric resonance. In these conditions 
a nonlinear coupling between the transverse and longi- 
tudinal modes plays important role. If the parametric 
growth of the amplitude of longitudinal modes advances 
their relaxation, after a noticeable reduction of b± one 
can expect a reverse nonlinear pumping of the transverse 
mode. As a result, the growth of b± starts and only 
later it will be succeeded again by the normal damping 
of the breathing mode. A specific feature of the paramet- 
ric resonance is its dependence on a product of the 
transverse and longitudinal modes. For the smaller mag- 
nitude of b± (t = 0), if a typical time for the pumping of 
the longitudinal mode proves to be comparable with the 
relaxation time or larger, the damping regime changes 
significantly and the amplitude of the breathing mode 
decays gradually in time. The picture described for the 
nonlinear damping allows to explain the results observed 
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in g(J. 

For the quantitative description of the picture of non- 
linear damping, we consider a simple model assuming 
that the parametric amplification in an elongated trap is 
experienced by a single longitudinal mode with frequency 
f2|| close to a half of the transverse mode frequency fl± 
and inverse relaxation time 7. We assume that the sys- 
tem is in the Thomas-Fermi regime and only long wave- 
length longitudinal sound modes are involved into the 
process. 

The Hamiltonian of interacting excitations in the Bose- 
Einstein condensate within second quantization for the 
representation of noninteracting phonons can be written 
as (see, @) 



H = H + H', 
H = J2 ?>Ms (c+c s + |) , H' = ±mJ d 3 xVh'V. 



(1) 



The operators of the alternating fraction of density h' 
and velocity V have the form 



n = ^[-2^4) (P.X.-C+X.), 

^ = E(^) 1/2 (^ s + a^). (3) 



(2) 



In these expressions noo is the magnitude of density 
no (x) at x = 0, Cq = fi/m, and m is the atom mass. The 
real eigenfunctions \s (%) for noninteracting phonons in 
a trap are solutions of the equation (cp.|Q) 



n2x.(^) + cgv (/ (x)V Xs (x)) = 0, 



(4) 



where / (x) = n (x) /n ao . 

Let us choose the breathing mode and one longitudi- 
nal mode in the excitation spectrum and, correspond- 
ingly, eigenfunctions Xj_(x) and X\\(%)- Accordingly, in 
this case there are two types of the vertices with in- 
tegrals like J = /d 3 xxjJx)Vx||(x)Vx||(x) and J' = 
J d 3 xx|| (x)Vxy (x)Vxj_(x) In the Thomas-Fermi approx- 
imation the function /(x) in Eq. ([J) equals 



m = 1 - w 



(5) 



where R± = y/2^/mu; 2 L and R\\ = ^J2fi/muj 2 are the 
transverse and longitudinal Thomas-Fermi radii. 

Let wave functions x s be normalized by the Thomas- 
Fermi volume V = ^TrR x R\\. Then, for function x± we 
find readily from Eq.(T 



X_L 




R 2 



(6) 



In first approximation in the limit uj\ ^> u| function X|| 
depends on z alone (e.g.,jllj). Hence, calculating integral 
J with use of @, we arrive at \J\ ~ uj\/\/Vcq. On the 
other hand, the estimate of J' yields \J'\ ~ lj?,/VVcq. 
In the strongly elongated trap we can retain only the 
terms with J in H'. As a result, in Hamiltonian H' 
there are only two terms from eight ones that can realize 
a quasi resonance coupling between these selected modes. 
Retaining only these terms, we represent the expression 
for H' in the explicit form 



H = ihgj(c+c+c ± - c+c^), 
1/2 

1 / hcf{il± \ 

9 = 




(7) 
(8) 



[cs,H]. 



In the Heisenberg representation 

dt Cs ~ ~Ji 

Using Eqs. ([[J) and (0), we find 

JjC|| = — ifi||C|| + 2gJctc x 
-iQ±c±_ - gJc..c.. 



(9) 



Let us go over into the classical Bose field in equations 
(^) . This is a conventional transformation in treating ex- 
citations in the Bose-Einstein condensate, see, e.g., po[| . 
In addition, we note that even for the extremely low tem- 
perature 40 nK p6[ , there are 5 phonons at the longitu- 
dinal mode and hundreds at the excited transverse mode. 
Within the framework of the standard procedure 

c,| -> b ]l e- i ^ t - v \0,c ± -» b±e- i ( Q±t - p ± ) . (10) 

Here b s and (p s are the slowly varying real quantities. 
Performing this transformation and separating the real 
and imaginary parts, we find 

d - - - - - 

— 6|| = -a&_L&||Cos</>-7(&|| - b||(0)), (11) 

di b± = 2< 
d hi 

at' 



-b± = —ab 2 cos <p, 



(12) 



-w = A^ + 26_l(1~ -=£-)sin<p. (13) 
ib ± 

Here we have singled out the initial value 6_i_(0) for 
breathing mode amplitude and introduced a relative am- 
plitudes b s = b s /b±(0) and the notation 

a = 2g\J\b±(p), (J<0). (14) 

The phase ip is defined as 

tp(t) = ip x (t) - 2y>||(t) + Am, Ail = tt± - 2Cl\\, (15) 



In the equation for the longitudinal ( p.l| ) component we 
have introduced the term taking into account effective 
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relaxation. Quantity b\\ (0) is the initial equilibrium value 



A set of Eqs. (11)-([L3|) describes nonlinear damping of 
the breathing mode in the general case. Let us consider 
first the initial period provided 6||(0) <C b±(0). Omitting 
the ratio (b\\/2b±) 2 and neglecting variation of b± in the 
Eq. ( |I3| ) we arrive to the equation with separable vari- 
ables. The solution of this equation can readily be found. 
For lAfil < 2a 



Afl tan % + 2a 
In 2 



Ail tan § + 2a + x 



xt + C, 



(16) 



where x — y (2a) — (Af2) . For xt ^> 1, phase <p 
tends to a constant value. For Ail — ► 0, ip — > 7r and 
in Eqs. (|ll|), ( |l^ ) cos<£> — > — 1. As a result, solution of 
Eq. ( p] ) can be represented as by oc exp[(a — 7)*]. An 
exponential growth with smaller decrement (|cos(/j| < 1) 
holds for finite Af2 if |AO| < 2a and a|cos</j| > 7. 
This is a typical parametric resonance resulting in the 
damping of the transverse oscillations in the condensate 



For an arbitrary time, the analysis requires a combined 
solution for all three nonlinear equations ( 11 )-(p~3|)- 

Let us find relation between the coefficients in this 
system and the physical parameters. The vibrational en- 
ergy, associated with the breathing mode, equals E v a w 
^fj,N(6R±/R±) 2 , where N is the total number of par- 
ticles. On the other hand, in accordance with defi- 
nition (fl"o|), one has E v u, = hil±b 2 L . Thus b±(0) — 
(4^N/3hQ ± ) 1/2 \SR±/R±\ , where \SR±/R±\ is an ini- 
tial value of a relative amplitude of the transverse oscil- 
lations. Estimating integral J at the vertex of the inter- 
action Hamiltonian (0) , we find J » 

Hence, determining coefficient a (]1< 
value of g (pf) , we have 



and involving 



a « 0.3cj 1 



SRi 



Ra 



(17) 



In derivation we used the following relations N — 
(2/5) noo^, 0|| Ktn±/2 = u±. 

The parametric resonance requires nonzero amplitude 
6||(0). For T = 0, b\\ (0) is the amplitude of zero-point 
oscillations. For T > Ml\\ we have 6|| (0) « (T/h^) 1 / 2 . 

Let us turn to the relaxation in the longitudinal sub- 
system. The familiar dispersion law for longitudinal ex- 
citations uj\ = w|fc(fc + 3)/4 [fl9"||,pl|| in combination with 
the mesoscopic character of the system results in a pecu- 
liar picture of relaxation. For a strongly elongated trap 
lo± ^> uju, frequency fiii is met with large value k ^> 1. 
This state can decay into numerous pairs of excitations 
with integer indices satisfying condition ki + &2 = k. It is 
interesting that the mismatch of resonances for different 
pairs changes only within interval Auj — u>h/2 ~ 3a>ii/4. 



Within the framework of the hydrodynamic approxima- 
tion these processes are described by the same Hamilto- 
nian H (|lj) involving (||) and (||) . Selecting a part of the 
interaction Hamiltonian responsible for the evolution of 
state k, we have 



C T, c fc 



ctc,_,c. 



(18) 



A direct estimate of vertex B k > yields 

By ~ 2i f (uJkUJk<u k -k>) 1/2 ■ 

In the usual conditions \B k >(k)\ <C Auj. At T = the 
excited /c-state spreads partially over k/2 pairs of states 
(for simplicity we imply that k is even) with the inverse 
time of reducing the density 



^ \BjAkf 
4^ Acj fc ' 



(19) 



All k states appeared spread, in turn, over a large 
number of states. In these conditions, taking also into 
account the scatter of Awy and vertices B k i, one can 
suppose that the probability for revival of fc-state with a 
noticeable amplitude is negligible. In essence, relaxation 
occurs for the time about r k (19). 

As concerns finite temperatures T > huik, the 
role of transitions with the involvement of low fre- 
quency states enhances. Hereupon, \B kl (k)\ 2 is multi- 
plied approximately by a factor T jhwy + T/fvjJk-k' — 
Tu> k / (hw k >uj k ^ k >) . Then 



1 



2fcw, 



(20) 



The parameter 7 in Eq. (|ll|) is w (2rfc) _1 . 

Let us consider a solution of the system of nonlinear 
equation s (|ll|) -(|r^). We choose the parameters close to 
that of (1J: ui = 1150 s -1 , oj\\ = 70 a" 1 , N = 10 5 , 
fi = 60 nK. The fmiteness of temperature T — 40 nK re- 
ally affects only the magnitude of parameter 7, see (po|). 
For these values of parameters k m 30 and coefficient 
7 lie within interval 5 10 s _1 . We take 7 = 7 s" 1 



for calculation. For initial condition \SR±/Rj 



0.09, 



coefficient a (|lj) takes value a w 30 s _1 . The spacing 
between the neighbor levels in the spectrum of longitu- 
dinal excitations is close to wm/2. Within the two-level 
model under consideration it is naturally to suppose that 
frequency wi, falling within this interval, is shifted to 
one of the levels. For calculating, we put AO, = 15 s _1 . 

In Fig.l the results of solving the system (pl])-([T^) 
are presented for the given values of a, 7 and Ail 
as solid lines. The time behavior of \SR±/R±\ = 
b± \SR±/R±\ is plotted in Fig. la. A drastic reduction 
of the breathing mode amplitude for the relatively short 



4 



0,10 



5a 

H 
*! 

■O 0,04 



: \ 


(a) 


: L 














(b) 







0,0 0,2 0,4 0,6 0,8 1,0 1,2 

time(s) 

Figure 1: (a) The amplitude of breathing mode \SR±/R± \ and 
(b) the amplitude of longitudinal mode b\\ as a function of time 
for two initial values \8R±/R± | = 0.09 (— ) and \SR±/R±\ = 
0.035 (---) 



order to make comparison more adequate. At the same a 
detailed quantitative comparison requires more elaborate 
model. 

It is worthwhile to note that an additional damping 
due to extrinsic reasons can be significant for the slow 
dynamics at small value of amplitude \8R±/ R±\ . How- 
ever, this damping will not play any noticeable role in the 
case of large magnitude of the amplitude. In this aspect 
the nonlinear damping is a purely intrinsic phenomenon. 

It seems that the results obtained have a relatively gen- 
eral character. In fact, the typical feature of all systems 
used for studying Bose-Einstein condensation in ultra- 
cold gases is their mesoscopic character and separation 
from environment. In these conditions the problem of 
the damping of coherent oscillations in the condensate at 
T = or, at least, T < \i <C T c acquires a special sound- 
ing. On one hand, this is connected with the problem 
of revealing the decay channel. On the other hand, pro- 
vided such channel exists, the energy is transferred to the 
discrete degrees of freedom. In this case the analysis of 
energy relaxation requires involving nonlinear coupling 
between modes for the both dissipationless and dissipa- 
tive kinetics. The problem considered in the present work 
is an instructive example of this general problem. 

This work is supported by INTAS (project 2001- 
2344) , by Netherlands Organization for Science Research 
(NWO), and by the Russian Foundation for Basic Re- 
search. 



interval of time is a direct manifestation of paramet- 
ric resonance. This is especially clear from the burst- 
like growth of the amplitude for longitudinal oscillations 
= b\\(t)\6Rj_/R±\ , as is seen in Fig. lb. At some 
time moment the phase ip reaches value Sir/ 2 (ir/2 for 
ACi < 0) and the reverse energy transfer from longitudi- 
nal into the transverse mode starts. In place of reduction, 
a nontrivial growth of the breathing mode amplitude 
starts. The growth is over when the phase reaches value 
5ir/2 (— 7r/2) and the normal trend recovers. Note that 
for the chosen magnitude of initial amplitude |<5i?j_/-Rj_| 
or parameter a, entire qualitative picture of nonlinear 
damping holds for a wide variation of Af2. 

However, the picture of damping changes drastically if 
the initial amplitude is noticeably smaller. In Figs, la 
and lb the same dependences for the case \SR±/R±\ Q = 
0.035 are shown with the dashed lines, the other param- 
eters being unchanged. We see a normal monotonic be- 
havior for the amplitude reduction \SR±/ R±\. As is seen 
from Fig. lb, the damping occurs in the conditions of the 
relatively small-expressed parametric resonance. Here- 
with, phase cp changes within the limited boundaries. 

The picture of nonlinear damping, found within the 
framework of solving a relatively simple model, explains 
the results observed in pi]. We have chosen intention- 
ally the magnitude of parameters close to that of |Q in 
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